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due  to  a  time  dependent  mechanism. 
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INTRODUCTION 


Since  the  early  sixties  there  has  been  considerable  interest  in  silicon 

nitride  as  an  engineering  material.  Sage  and  Histed^  reviewed  some  of  the 

applications  for  silicon  nitride  and  suggested  the  use  of  silicon  nitride  in  gas 

turbines.  An  increase  of  100°C  in  turbine  stator  blade  temperature  was 

(2) 

forecast.  Glenny  and  Taylor  reviewed  the  material  properties  of  silicon 

nitrides  and  noted  the  high  strength  at  elevated  temperatures  of  silicon  nitride. 

(3) 

Evans  and  Wiederhorn  observed  the  phenomenon  of  static  fatigue  in  hot 
pressed  silicon  nitride,  and  noted  that  this  presented  a  design  problem  in  the 
application  of  the  material  in  highly  stressed  applications  at  elevated 
temperatures. 

The  phenomenon  of  static  fatigue  is  well  known  in  metals  where  it  may  be 
caused  by  slow  crack  growth.  Siverns  and  Price  ^  observed  that  the  crack 
growth  rate  in  steel  at  elevated  temperatures  was  correlated  with  the  power 
law  relationship 

V  =  AKjn  .  (1) 

Evans  and  Johnson^  suggested  that  this  equation  could  be  applied  to  ceramics 
and  that  the  cause  of  slow  crack  growth  was  grain  boundary  sliding.  This  is  due 
to  the  glassy  phases  in  the  grain  boundary  softening  at  elevated  temperature 
and  allowing  the  crystalline  grains  of  silicon  nitride  to  slide  relative  to  one 
another  at  the  tip  of  the  crack  where  the  stresses  are  the  highest.  The  sliding 
of  the  grains  causes  the  crack  to  grow  and  continued  application  of  stress  and 
temperature  could  cause  catastrophic  failure^. 

Williams  and  Evans^  proposed  the  double  torsion  test  to  measure  the 
values  of  the  crack  velocity  exponent,  and  the  premultiplier  in  Eqn.  1.  The 
double  torsion  test  is  a  simple  test  in  concept  that  loads  a  pre-cracked  specimen 
in  a  manner  that  the  stress  intensity  factor  remains  constant  as  the  crack 
grows.  This  simplifies  the  analysis  and  makes  it  possible  to  determine  A  and  n 
from  the  load  versus  time  curve. 

(8) 

Charles  proposed  stress  rate  testing  to  determine  the  crack  velocity 
exponent  and  demonstrated  its  use  on  glass.  In  stress  rate  testing  MOR  bars  are 
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fractured  at  two  different  stress  rates  to  obtain  two  pairs  of  data  -  the  fast 
fracture  stress  and  its  corresponding  stress  rate.  From  the  pairs  of  data  the 
slow  crack  growth  parameters  are  calculated. 

Several  experimenters  have  conducted  stress  rupture  tests  on  ceramic 
(9) 

materials.  Trantina  conducted  stress  rupture  tests  on  Norton's  NC-132 
material,  a  hot  pressed  silicon  nitride.  His  tests  were  conducted  in  three-point 
bending  on  bars  of  0.9  inch  span  with  a  cross  section  of  0.1  inch  x  0.1  inch. 
Quinn^^  conducted  stress  rupture  tests  on  several  different  structural  ceramic 
materials.  His  tests  were  conducted  in  four-point  bending  on  a  fixture  with  a 

I. 5  inch  outer  span,  a  0.75  inch  inner  span,  and  a  bar  with  a  0.110  inch  base  and 
a  0.085  inch  height. 

Previous  experimenters  conducted  their  tests  on  small  parts  with  very 
small  volumes  compared  to  the  volume  of  typical  gas  turbine  components. 
Previous  tests  were  also  conducted  at  uniform  temperatures,  but  gas  turbine 
components  have  large  temperature  gradients.  This  suggested  the  need  for 
generating  a  time  to  failure  database  on  a  component  the  size  of  a  typical  gas 
turbine  component  with  temperature  and  stress  distributions  similar  to  gas 
turbine  components.  This  database  could  then  be  used  as  a  test  case  for 
evaluation  of  models  of  failure  and  computer  programs  used  to  predict  time  to 
failure. 

The  hot  spin  rig  shown  in  Fig.  1  has  accumulated  many  hours  of  testing 
ceramic  gas  turbine  rotors  at  elevated  temperatures  and  at  angular  speeds  equal 
to  gas  turbine  speeds.  The  rig  proved  very  reliable  under  these  conditions, 
therefore,  it  was  selected  to  conduct  the  tests  required.  Operation  and  details 
of  the  rig  have  been  previously  described^  A  rotating  disc  was  designed  to 
fail  due  to  time  dependent  mechanism  utilizing  available  techniques  and  models. 
The  rig  was  developed  to  conduct  the  tests,  and  an  initial  calibration  test  was 
conducted. 

II.  DESIGN  OF  DISK 

The  finite  element  grid  of  the  rotating  components  used  to  calculate  the 
temperature  distribution  in  the  spin  disc  is  shown  in  Fig.  2.  The  physical  parts 


-2- 


represented  by  the  finite  elements  are  shaded  and  labeled  for  identification. 
The  finite  elements  that  compose  the  model  are  rectangular  or  triangular  in 
shape.  The  material  thermal  properties  used  in  the  analysis  are  given  in  Table  I. 


The  temperatures  within  the  rotating  components  were  determined  by 
these  boundary  conditions  and  the  thermal  conductivities  of  parts  in  the 
assembly.  Heat  was  applied  to  the  front  face  and  rim  of  the  spin  disc  with 
propane  gas  burners.  The  temperature  of  the  front  face  was  measured  with  a 
radiation  pyrometer  and  the  measured  temperatures  were  used  as  input  to  the 
analysis.  The  tie  bolt  was  cooled  with  air  as  shown  on  Fig.  2.  The  cooling  air 
mass  flow  rate  and  associated  heat  transfer  coefficients  formed  the  boundary 
conditions.  A  portion  of  the  tie  bolt  cooling  air  was  directed  to  the  metal 
curvic  adaptor  for  cooling  and  the  mass  flow  rate  and  heat  transfer  coefficients 
formed  the  boundary  conditions  for  the  exterior  of  the  tie  bolt,  nut,  and  metal 
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curvic  adaptor.  The  downstream  side  of  the  spin  disc  and  the  exterior  of  the 
ceramic  curvic  spacer  were  heated  with  hot  gases  and  air  swirling  in  the  rig, 
thus  forming  another  boundary  condition.  The  labyrinth  and  bearing  region  on 
the  shaft  form  the  last  boundary  condition.  As  part  of  the  model,  radiation  heat 
transfer  was  included  between  the  spin  disc  bore,  ceramic  curvic  spacer  bore 
and  the  tie  bolt. 


GAS  FLOW 


1  SPIN  DISC 

2  CURVIC  * 

3  CERAMIC  CURVIC*SPACER 

4  SHAFT 


Figure  2-Rotating  components  heat  transfer  model. 

Calculated  temperature  distributions  were  compared  to  results  from 
temperature  sensitive  paints  and  adjustments  to  the  model  were  made  to  obtain 
agreement.  Temperature  measurements  of  the  air  in  the  cooling  streams  were 
also  used  as  input  to  the  model  and  as  a  check  on  the  validity  of  the  model;  for 
example  the  temperature  of  the  tie  bolt  cooling  air  was  measured  as  it  entered 
the  labyrinth  and  as  it  left  the  tie  bolt.  The  temperature  of  the  air  in  the 
labyrinth  region  was  used  as  input  to  the  model  and  the  air  exit  temperature 
was  compared  to  the  calculated  value  as  a  check  on  the  model. 
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Figure  3  illustrates  the  finite  element  grid  used  to  calculate  the  stresses 
in  the  spin  disc.  The  temperature  distribution  was  derived  from  the  heat 
transfer  model  and  transferred  to  the  stress  finite  element  grid.  The 
temperature  contour  plot  at  the  test  design  conditions  is  illustrated  in  Fig.  4. 
The  combined  thermal  and  centrifugal  stresses,  at  this  test  condition,  are 
illustrated  in  Figs.  5  through  9.  These  stress  plots  show  the  maximum  principal 
tensile  stress,  the  tangential  stress,  the  radial  stress,  the  axial  stress,  and  the 
maximum  shear  stress.  The  elastic  properties  used  in  the  stress  calculations  are 
given  in  Table  11. 


Figure  3-  Disc  finite  element  grid. 
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Figure  8  -  Axial  stresses  [psi]  in  disc  combining  thermal  and  centrifugal 
stresses  at  50,000  rpm. 


Figure  9  -  Maximum  shear  stresses  [psi]  in  disc  combining  thermal  and 
centrifugal  stresses  at  50,000  rpm. 
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(12) 

The  reliability  versus  time  formula  used  in  this  report  is  shown  in  Eqn. 

(1).  Definition  of  the  symbols  is  given  in  appendix  B. 

n-2  —i  m 

m  a  nt  n-2 


R  =  exp 


On  b 

Rff 


°<P~ 


'B 


(1) 


During  this  work,  Eqn.  1  was  reduced  to  a  simpler  form  shown  as  Eqn.  2.  The 
detailed  algebraic  steps  are  given  in  Appendix  C. 


m 


R  =-R 


/.  (J  ti  n-2 
(1+  -=-) 


it 


B 


(2) 


Equation  2  is  much  simpler  for  computation  and  was  programmed  for  digital 
computer  use  to  calculate  time-dependent  reliabilities  for  steady  state 
temperature  and  leading  conditions.  The  equation  was  evaluated  for  each 
element  of  the  finite  element  grid  of  Fig.  3  based  on  the  maximum  principal 
stress,  the  Weibul!  modulus,  the  crack  velocity  exponent,  the  value  of  B,  and  the 
fast  fracture  reliability  of  the  element.  The  Weibuli  modulus,  crack  velocity 
exponent,  and  the  value  of  B  were  linearly  interpolated  from  a  table  of  these 
parameters  versus  temperatures,  based  on  the  centroidal  temperature  of  the 
element. 


Equation  2  shows  that  when  time  equals  zero  the  time  dependent 

reliability  is  simply  the  fast  fracture  reliability.  The  fast  fracture  reliability 

was  evaluated  by  integrating  the  normal  stress  around  the  unit  sphere  in  each 

(13  14) 

finite  element  and  evaluating  the  Weibuli  integral  for  each  element  *  . 

Two  data  bases,  double  torsion  data  and  stress  rate  data,  were  available 
for  the  calculation  of  time  dependent  reliabilities.  Both  data  bases  yield  the 
two  parameters  "B"  and  "n"  essential  for  the  evaluation  of  Eqn.  2.  Figure  10 
shows  the  calculation  flow  path  utilizing  double  torsion  data.  Fast  fracture 
testing  produced  a  collection  of  fast  fracture  data  which  was  analyzed  using 
maximum  likelihood  estimation  methods1^"5)  to  determine  the  characteristic 
modulus  of  rupture  and  the  Weibuli  modulus.  The  stress  distributions  along  with 
the  Weibuli  parameters  were  input  to  a  computer  program  which  evaluated 
Weibull's  formula  for  fast  fracture  reliability,  which  in  turn  was  input  to 
another  computer  program  which  calculated  time  dependent  reliability.  From 


double  torsion  testing  the  crack  velocity  exponent,  the  premultiplier  W,1  and  the 

critical  stress  intensity  factor  were  obtained.  With  these  parameters  and  the 

value  of  the  stress  intensity  factor  coefficient,  "Y",  the  value  of  "B"  was 
(12) 

calculated  .  The  value  of  "B"  and  crack  velocity  exponent  "n"  were  input  to 
a  computer  program  that  evaluated  the  reliability  versus  time. 


FAST  FRACTURE 
TESTING 


DOUBLE  TORSION 
TESTING 


n,A  Kjp 


RELIABILITY  VS  TIME 


Figure  10  -  Calculation  flow  path  utilizing  double  torsion  data 


Figure  11  shows  the  calculation  flow  path  utilizing  stress  rate  testing  data 
to  obtain  reliability  versus  time.  Stress  rate  testing  determined  pairs  of  fast 
fracture  strength  and  the  corresponding  value  of  stress  rate.  The  value  of  "B" 
and  the  crack  velocity  exponent  "n"  were  calculated  from  the  pairs  of  strength 
and  stress  rate  with  the  formulas  shown  on  the  figure  The  fast  fracture 
calculation  was  the  same  as  that  used  with  double  torsion  data. 

The  calculation  procedures  shown  in  Figs.  10  and  1 1  were  used  to  design 
the  disc.  The  design  was  arrived  at  by  checking  trial  designs  with  the  reliability 
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versus  time  calculation  procedures  until  a  satisfactory  design  was  obtained. 
The  relibability  versus  time  results  are  shown  on  Fig.  12.  Note  that  the  two 
data  bases,  double  torsion  and  stress  rate,  give  widely  different  results. 


FAST  FRACTURE 
TESTING 


STRESS  RATE 
TESTING 

®i,0ii  i®2,<ri2 


RELIABILITY  VS  TIME 


Figure  11  -  Calculation  flow  path  utilizing  stress  rate  data 
III.  DISC  MACHINING 

Life  prediction  discs  were  machined  from  six  inch  square,  one  and  three- 
eighths  inch  thick  billets  of  hot  pressed  silicon  nitride  (NC132).*  The  four  inch 
diameter  discs  were  made  in  several  steps  beginning  with  a  slicing  operation  in 
which  an  octagonal  block  was  cut  out  of  each  billet.  This  piece  was  ground  to 
the  approximate  thickness  of  the  finished  disc  and  a  center  hole  or  bore  was 
then  core  drilled  through  with  a  diamond  core  drill.  The  bore  was  then  ground 


*  Certified  Billets  from  Norton  Company 


-  11  - 


to  finish  size  allowing  the  disc  to  be  centrally  mounted  for  further  grinding 
operations.  Most  of  the  grinding  was  done  on  a  Brown  and  Sharp  Universal 
grinder.  Initial  material  removal  was  accomplished  using  120-150  grit  diamond 
plate  and  a  commercial  water  soluble  coolant  (X50-2  Chem-trend).  Final 
dimensions  were  produced  with  200-250  grit  wheels. 


TIME  (HOURS) 

Figure  12  -  Calculated  reliability  versus  time 


The  disc  face  contours  were  produced  by  specially  shaped  contour  grinding 
wheels  shown  in  Fig.  13.  These  wheels  were  rotated  at  12-13000  rpm  and  were 
fed  into  the  disc  side  as  shown  in  Fig.  14  in  a  three  step  operation.  The  disc  was 
rotated  at  100  rpm  in  the  opposite  direction  to  the  grinding  wheel  rotation. 
This  allowed  control  of  the  grinding  and  stability  of  the  set-up. 
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Figure  13  -  Contour  grinding  wheel 


3  STEP  PROCEDURE  FOR  DISC  CONTOURING 


SURFACE 


SPECIFIED 

DIMENSION 


GRINDING  WHEEL  BROUGHT 
AXIALLY  INTO  CONTACT  WITH 
SURFACE  V  AND  RADIALLY 
INTO  BASIC  HUB  DIAMETER. 


GRINDING  WHEEL  ADVANCED 
PARALLEL  TO  DISC  AXIS. 


FINISHED 


STEP  3 _ 

GRINDING  WHEEL  ADVANCED 
PERPENDICULAR  TO  GRINDING 
WHEEL  AXIS  TO  COMPLETE 
CONTOURING. 


MATERIAL  IN  CROSS  HATCHED 
VOLUME  REMOVED. 


Figure  14  -  Procedure  for  disc  contouring 
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All  grinding  tools,  including  contour  wheels,  were  diamond  coated  and 
final  polishing  of  the  disc  contours  was  done  with  diamond  pastes  on  a  brass 
lap.  The  lap  was  in  the  form  of  a  ball  on  a  1/4"  diameter  shank  and  the  lap  was 
driven  by  a  hand  held  electric  drill.  The  disc  was  slowly  rotated  on  an  arbor  and 
the  contour  was  polished  to  a  1-2/xinch  finish  by  successively  using  diamond 
pastes  of  30-40,  6-12  and  1-5  micron  grit. 

Finally  the  curvic  couplings  used  to  mount  the  disc  to  the  shaft  were 
ground  by  diamond  coated  grinding  wheels  on  special  purpose  curvic  grinding 
equipment.  The  discs  were  then  zygloed  for  possible  cracks  caused  by  grinding 
and  then  dimensionally  inspected.  Figure  15  shows  the  completed  disc  ready  for 
testing. 


Figure  15  -  Disc  ready  for  testing 
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IV.  RIG  DEVELOPMENT 


The  hot  spin  rig  was  originally  developed  to  conduct  tests  on  bladed  discs. 
Considerable  experience  had  been  obtained  on  bladed  discs  at  a  platform  tern- 
perature  of  1800°FU  .  In  addition  25  hours  of  testing  had  been  conducted  at  a 
platform  temperature  of  2000°F  at  a  speed  of  50,000  rpm^^.  Three  hours  of 
testing  had  been  conducted  at  a  platform  temperature  of  2200°F  and  speeds 
above  30,000  rpm.  Since  no  direct  experience  on  unbladed  discs  at  the 
temperatures  required  for  the  Hot  Spin  Disc  Life  Program  existed,  a 
development  program  was  necessary. 

Tests  on  actual  Life  Prediction  discs  began  in  mid  July  of  1980,  with  the 
initial  attempts  to  run  disc  NC1316  at  2300°F  rim  temperature.  The  original 
burner  configuration  with  5  straight  burners  was  found  to  give  a  maximum  rim 
temperature  of  2130°F.  A  calibrated  radiation  pyrometer,  Ircon  model  300T5, 
was  sighted  through  a  hole  in  the  rig  wall  to  determine  the  actual  disc 
temperature.  A  slight  temperature  dependence  on  speed  was  found  but  it  was 
obvious  that  the  original  system  would  not  produce  the  required  2300°F  rirn  on 
an  unbladed  disc. 

The  first  modification  to  the  burner  system  was  to  increase  the  number  of 
burners  from  5  to  12.  Tests  with  this  set-up  produced  maximum  rim 
temperatures  of  2170°F. 

Tests  with  individual  burners  showed  that  the  gas  flame  temperature  was 
running  at  3100  to  3200°F;  however,  the  burners  were  located  such  that  only 
part  of  the  gas  flame  touched  the  disc  rim.  The  flames  did  not  impinge  directly 
onto  the  disc  rim.  It  was  decided  to  modify  two  burners  so  that  the  flame  would 
impinge  directly  on  a  portion  of  the  rim  of  the  front  face  of  the  disc.  These 
burners  were  cut  and  rewelded  at  an  angle  of  10°  and  the  total  number  of 
burners  was  reduced  to  five.  In  addition,  as  shown  in  Fig.  16,  the  nose  cone 
portion  of  the  rig  was  redesigned  and  fabricated.  The  original  design  ducted  the 
bolt  cooling  air,  exhausted  through  the  end  of  the  bolt,  back  over  the  curvic 
adaptor  and  up  the  forward  face  of  the  disc.  This  600°F  air  produced  a  veiling 
as  well  as  a  cooling  effect.  The  redesigned  parts  ducted  this  cool  exhaust  air 
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out  of  the  rig  through  a  hole  along  the  center  of  the  nose  cone  exhausting  to 
ambient.  Also  some  of  the  rig  insulation  behind  the  disc,  which  ran  at  very  high 
temperatures,  was  increased  in  area  to  increase  the  hot  surfaces  radiating  heat  ' 

to  the  rear  face  of  the  disc. 


EXHAUST  TO  REGULATION 


Figure  16  -  Tie  bolt  cooling  air  distribution 

This  design  was  tested  in  mid  October  of  19S0.  The  first  tests  were  run 
without  restricting  the  bolt  exhaust  air.  A  rim  temperature  of  2370°F  was 
easily  obtained;  however,  the  jet  effect  of  the  bolt  exhaust  air  drew  hot  gases 
down  the  disc  face  into  the  curvic  area  and  overheated  the  curvic  adaptor.  This 
over-heating  softened  the  adaptor  and  caused  the  rotor,  NC1316,  to  fail  at 
24099  rpm.  Apparatus  was  added  externally  to  the  rig  to  allow'  modulation  of 
the  bolt  exhaust  air.  This  caused  some  of  the  bolt  cooling  air  to  flow  over  the 
metal  curvic  adaptor  maintaining  it  at  an  acceptable  temperature.  A 
thermocouple  was  installed  in  the  adaptor  cooling  air  flow  stream  to  monitor 
the  temperature. 
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These  modifications  completed  the  rig  development,  the  details  of  the  tie 
bolt  cooling  air  distribution  are  shown  in  Fig.  16.  This  illustrates  how  the  tie 
bolt  cooling  air  is  split  so  that  a  portion  of  the  air  flows  over  the  metal  curvic 
adaptor  for  cooling.  The  overall  schematic  of  the  final  rig  configuration  is 
shown  in  Fig.  17.  The  modifications  to  the  burner  tubes,  nose  cone,  and 
insulation  may  be  seen  by  comparing  Fig.  17  to  Fig.  1. 


CONSTANT  AM 
PLOW  PLENUM 


CEAAMK 

SPACCT 


TUT  DISC 


SUB- ATMOSPHERIC 
EXHAUST 


FUSED  SILICA  SHROUD 


Figure  17  -  Final  rig  configuration 


V.  INITIAL  CALIBRATION  TEST 

An  initial  calibration  test  was  conducted  on  disc  NC  1318  to  verify  the 
analytical  calculations  and  to  determine  if  failure  would  occur  within  a 
reasonable  time  frame.  The  rig  configuration  used  is  shown  in  Fig.  17.  The 
steady  state  temperatures  selected  were  the  same  as  those  shown  in  Fig.  4  with 
a  corresponding  speed  of  50,000  rpm.  The  steady  state  conditions  were  arrived 
at  by  igniting  the  burners  with  the  disc  rotating  at  2000  rpm.  The  disc  was  then 
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accelerated  to  10,000  rpm.  The  speed  and  temperature  were  stabilized  at 
10,000  rpm.  The  disc  was  then  accelerated  to  24000  rpm,  where  the  speed  and 
temperature  were  stabilized.  The  temperature  gradient  was  established  at 
24000  rpm  so  that  the  thermal  stress  distribution  was  established  and  the 
relative  motion  between  rotating  components  would  be  minimizea  during  the 
acceleration  to  50,000  rpm.  After  the  temperature  gradient  stabilized,  the  disc 
was  slowly  accelerated  to  50,000  rpm.  At  50,000  rpm  the  temperatures  and  air 
flows  were  stabilized  and  they  remained  constant  for  the  rest  of  the  test. 

Failure  occurred  eight  and  one-half  hours  after  reaching  50,000  rpm.  No 
indications  of  problems  occurred  on  any  of  the  instrumentation  used  to  monitor 
the  test,  which  indicates  that  the  test  rig  was  performing  properly.  After  the 
test,  the  rig  was  disassembled  and  the  parts  carefully  examined  to  determine  if 
failure  occurred  due  to  slow  crack  growth.  The  stress  analysis  indicated  that 
the  most  severe  combination  of  temperatures  and  stresses  was  in  the  neck 
region.  Figure  18  is  a  photograph  of  the  hub  region  after  the  test,  showing  that 
the  hub  remained  intact.  This  shows  that  the  failure  did  not  initiate  in  the  hub 
region.  Figure  19  shows  a  piece  of  the  neck  region  showing  the  fracture  and  a 
crack  that  developed  parallel  to  the  main  fracture.  Figure  20  shows  this  piece 
rotated  so  that  the  main  fracture  surface  can  be  seen.  It  is  believed  that  the 
pattern  of  this  crack  indicates  slow  crack  growth  occurring  normal  to  the  radial 
stress  in  the  neck  region.  The  conclusion  was  that  the  test  procedure  and 
conditions  were  developed  to  the  point  where  discs  could  be  tested  to  failure 
under  conditions  of  slow  crack  growth. 


Figure  18  -  Hub  region  of  the  disc  after  the  test 
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Figure  19  -  View  of  the  crack  in  the  neck,  normal  to  the  machined  surface 


Figure  20  -  View  of  the  crack  in  the  neck  showing  adjacent  fracture  surface 


VI.  SUMMARY 


A  disc  was  designed  to  fail  under  the  conditions  of  slow  crack  growth. 
Several  discs  were  fabricated.  The  existing  hot  spin  rig  was  modified  and 
developed  to  test  the  disc  under  the  design  conditions.  An  initial  calibration 
test  was  successfully  conducted  to  verify  the  design. 

The  literature  was  searched  for  slow  crack  growth  material  properties 
required  to  calculate  reliability  versus  time. 

Reliability  versus  time  was  calculated  for  the  disc  using  stress  rate  and 
double  torsion  data. 

VII.  RECOMMENDATIONS 

A  series  of  discs  should  be  tested  under  the  conditions  documented  in  this 
report  to  experimentally  determine  reliability  versus  time.  The  experimental 
data  so  determined  could  then  be  compared  to  the  analytical  results  presented 
here,  this  comparison  would  determine  which  data  basegives  the  most  accurate 
reliability  predictions.  The  experimental  data  would  also  be  useful  in  verifying 
mathematical  models  of  time  dependent  failure  of  structural  ceramics  and 
associated  computer  programs  incorporating  these  models. 
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Appendix  A 


Calculation  of  Constant  "B" 

The  calculation  of  time  dependent  reliabilities  is  based  on  the  power 
relation  of  crack  velocity  as  a  function  of  stress  intensity  factor. 

V  =  AK"  A-l 

The  values  of  A  and  n  are  determined  experimentally,  so  a  literature  search  was 
conducted  to  find  published  values  to  utilize  for  design  and  analysis.  Table  III 
gives  the  experimental  values  that  were  found.  The  experimental  values  were 
published  in  graphical  or  numerical  form  and  the  manner  of  presentation  is 
shown  on  the  table.  In  the  cases  of  graphical  publication  the  numerical  values 
were  found  by  taking  values  of  velocity,  "V",  and  stress  intensity  factor,  "Kj", 
from  the  graph  and  substituting  into  Eqn.  A-l.  It  should  be  noted  that 
considerable  error  in  the  determination  of  A  and  n  may  result  from  this 
procedure.  Cargill  and  Annis  published  their  results  graphically;  however,  they 
furnished  the  values  of  A  and  n  determined  by  least  squares  fitting  of  their 
data.  The  A  and  n  values  of  Annis  and  Cargill^  ^  and  Govila^1^  were  used  in 
this  program. 

The  values  of  A  and  n  available  are  published  for  930°C,  1200°C,  and 
1300°C.  Reliability  calculation  requires  values  for  A  and  n  at  the  centroidal 
temperature  of  each  finite  element  of  the  finite  element  grid  shown  in  Fig.  2. 
To  obtain  the  required  values,  linear  interpolation  was  used. 

The  critical  stress  intensity  factor  for  the  material  is  an  important 

parameter  when  double  torsion  data  is  used  to  calculate  time  dependent  relia- 
21) 

bility.  Larsen  presented  graphical  data  for  hot  pressed  silicon  nitride  and 
the  numerical  values  used  are  shown  in  Table  IV.  Also  required  for  the  use  of 
double  torsion  data  is  an  assumption  about  the  shape  of  the  flaw,  so  that  the 
stress  intensity  factor  coefficient  may  be  determined.  For  this  program  the 
flaws  were  assumed  to  be  penny  shaped  internal  flow  with  a  stress  intensity 
factor  coefficient^^  of  2/^  . 


As  shown  in  Fig.  9  the  value  of  "B"  must  be  computed  before  reliability 
versus  time  can  be  calculated  and  the  following  equations  show  the  evaluation 
of  "B"  at  1200°C  with  the  required  material  and  geometric  constants  and  their 
sources. 


n 

=  27.4 

Reference 

19 

A 

=  1.58xl0“23  — 

i 

19 

KIC 

sec 

=  7.2  MPa /In 

(MPa  /m)n 

21 

Y 

=  2//r 

22 

(n-2)Y2  AKIcn‘2 

B  = _ 2  _ _ 

(27.4-2)  (-  )2  1.58  x  10‘23  x  (7. 2)27A~2 

sa 

B  =  0.000182  MPa2  -  hr. 

An  alternate  data  base  to  use  for  predicting  time  dependent  failure  is 

(23) 

stress  rate  data.  The  NC-132  stress  rate  data  is  given  in  Table  V 
Reference  23  gives  data  for  several  stress  rates,  but  the  lowest  available  stress 
rate  and  corresponding  fracture  strength  was  chosen  on  the  assumption  that 
slow  stress  rate  result  was  most  applicable  for  predicting  stress  rupture  results. 
As  shown  on  Fig.  10  the  value  of  "B"  must  be  evaluated  before  reliability  versus 
time  can  be  calculated  and  the  following  equations  show  a  typical  calculation. 


Temperature  =  1204°C 

n  =  19.9 

cr{  =  485  MPa 

a  =  0.8  MPa/min. 

af3 

B  =  a(n+l) 

a  _  (485)3  hr. 

D  “  0.8  x  (19.9+0  60  min. 

B  =  113720  MPa2  -  hr. 


Reference 

(23) 

(23) 

(23) 
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APPENDIX  B  NOMENCLATURE 


A  Premultiplier  in  crack  velocity  equation,  meter/second  (MPa-m)n). 

2 

B  Constant  in  reliability  versus  time  equation,  MPa  -hr. 
i  ith  value  of  test  in  a  set  of  tests. 

Kj  Stress  intensity  MPa  /m. 

Kj£  Critical  stress  intensity,  a  material  parameter,  MPa  /rn. 
n  Crack  velocity  exponent. 

N  Number  of  tests  in  a  set. 

m  Weibuli  modulus. 

MORq  Characteristic  modulus  of  rupture,  MPa 
R  Reliability  as  a  function  of  time 

Rjj  Fast  fracture  reliability 

a  Stress,  MPa. 

Og  Characteristic  strength  of  the  structure,  MPa. 

Oj  Fracture  stress,  MPa. 

a  Stress  rate,  MPa/min. 

j0  Weibuli  parameter,  MPa/(m^)m. 

T  Temperature,  °C. 

t  Time,  hours. 

V  Volume,  m^. 

V  Crack  velocity,  m/sec. 

V  Stress  intensity  factor  coefficient,  non-dimensional. 
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APPENDIX  C 

SIMPLIFICATION  OF  RELIABILITY  VERSUS  TIME  FORMULA 


Reference  12  gives  equation  1  for  calculating  reliability  versus  time  for  an 
element  of  a  structure. 


This  appendix  shows  the  algebra  required  to  simplify  Equation  1  to  Equation  2. 

m 


The  nomenclature  is  given  in  Appendix  B. 

First  use  Equation  3  which  is  the  fast  fracture  reliability  of  any  element. 


This  can  be  rearranged  to  give  Equation  4. 


Now  one  group  of  terms  in  Equation  1  can  be  written  as'  Equation  5« 


Substitute  Equation  4  into  Equation  5  to  yield  Equation  6. 
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Appendix  c  continued 


I 

l 

Sk.. 


How  use  Equation  6  to  write  liquation  I  as  Equation  7» 


(7) 


(8) 


Using  the  laws  of  exponent^  Equation  8  become st 


Using  the  laws  of  logarithms, Equation  9  becomes: 


From  the  laws  of  exponents,  using  a5^  =  (aX)  ,  and  from  the  laws  of 
logarithms  e^3  =  b,  Equation  10  becomes: 


Equation  11  simplifys  the  computation  of  reliabilities  considerably  from 
the  use  of  Equation  1.  Also,  it  simplifies  the  understanding  of  how  each 
parameter  R^,  o,  B,  t,  m,  and  n  affect  the  reliabilities  calculated  from 
Equation  1. 
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TABLE  I 


Thermal  Properties 


Temperature 

Thermal 

Specific 

°F 

Cor^ptivity 

hr  ft  °F 

Heat 

Btu/lb  -  F 
m 

Spin  Disc  and  Ceramic 

70. 

17. 

.18 

Curvic  Spacer  Material 

500. 

15. 

.23 

Hot  Pressed  Silicon  Niyide 
Density  =  3.18  gm/cm"5 

1000. 

13. 

.26 

1500. 

11. 

.29 

2000. 

9.2 

.33 

2500. 

8.0 

.32 

Nut  Material 

0. 

6.0 

.10 

Inconel  X-750  , 

Density  =  8.30  gm/crn 

3000. 

20.5 

.18 

Tie  Bolt  Material 

0. 

6.0 

.10 

Inconel  X-718  3 

Density  =  8.19  gm/crri 

3000. 

21.8 

.10 

Compression  Member 

0. 

6.5 

.11 

Material 

A-286  , 

Density  =  7.91  gm/cnri 

3000. 

25.2 

.11 

Shaft  Material 

0. 

21.7 

.11 

AMS  6265  3 

Density  =  7.83  gm/cnri 

3000. 

21.7 

.11 

Curvic  Material 

70. 

3.99 

.18 

500. 

3.87 

.23 

Density  =  3.18  gm/cm 

1000. 

3.72 

.26 

1500. 

3.54 

.29 

2000. 

3.33 

.33 

2500. 

3.16 

.32 

Air  Temperature 

Thermal 

Conductivity 

Specific 

Heat 

Density 

Btu 

Btu 

Ibm 

hr  ft  °F 

Ibm 

T 

it* 

0 

.013 

.239 

.0863 

100 

.015 

.240 

.0680 

500 

.023 

.247 

.0392 

800 

.028 

.256 

.0299 

1000 

.032 

.262 

.0258 

1500 

.040 

.276 

.0192 

2000 

.047 

.286 

.0153 

2500 

.051 

.292 

.0127 

3000 

.054 

.297 

.0108 
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TABLE  II 


Elastic  Properties 


NC-132 

Temperature 

Young's 

Poisson's 

Shear 

Coefficient 

Hot  Pressed 

Silicon  Nitride 

Modulus 

Ratio 

Modulus 

of  Thermal 

Expansion 

in/in/°F 

Density  =  3.18  gm/cm3 

°F 

xlO^psi 

xlO^psi 

xIO"6 

78 

41.9 

.219 

17.1 

.74 

500 

41.9 

.203 

17.5 

.97 

1000 

41.8 

.190 

17.5 

1.32 

1500 

40.6 

.189 

17.1 

1.53 

2000 

38.7 

.194 

16.1 

1.68 

2500 

35.4 

.194 

14.3 

1.82 

I 
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TABLE  III 


. .  a-s- *r *t" ■  ■  nH* -'v-~ tn*— --  «••--  •  --v.  •'#■■  <•-•  <■ 


Crack  Velocity  Parameters 


Material 

Temper¬ 

ature 

°C 

HPSN 

NC-132 

930 

HPSN 

NC-132 

1200 

HPSN 

NC-132 

1300 

HPSN 

NC-132 

1350 

HPSN 

NC-132 

1400 

HPSN 

NC-132 

1200 

HPSN 

NC-132 

1400 

HPSN 

HS-130 

1250 

HPSN 

HS-130 

1300 

HPSN 

HS-130 

1350 

HPSN 

HS-130 

1400 

HPSN 

HS-130 

1200 

HPSN 

HS-130 

1400 

HPSN 

W 

1400 

Crack 
Velocity 
Exponent 
A  n 


6.18xl0"56 

74.7 

1.58xl0"23 

27.4 

4.3xl0~9 

5.6 

7.9xl0-9 

5.6 

lxlO”7 

5.25 

1.21xl0-31 

35.4 

9x1  O'8 

5.8 

1.41xl0~14 

10.5 

9.31X10"14 

11.5 

9.83xl0-13 

10.2 

3.05xl0"12 

10.6 

2.56xl0~59 

73.3 

8.22xl0'8 

4.93 

5xl0"14 

6 

Form  of 
Data 

Presentation 

Reference 


19 

Graphical 

19 

Graphical 

20 

Numerical 

20 

Numerical 

20 

Numerical 

24 

Graphical 

24 

Graphical 

25 

Graphical 

25 

Graphical 

25 

Graphical 

25 

Graphical 

25 

Graphical 

25 

Graphical 

26 

Numerical 
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TABLE  IV 


Critical  Stress  Intensity  Factor 

Temperature 

Stress  Intensity 

Factor 

°C 

MPa  ^7n 

18 

5. 

930 

5. 

1200 

7.2 

1300 

11.0 

Reference  21 
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TABLE  V 


STRESS  RATE  DATA  OF 
NORTON  NC-132  HPSN 
Characteristic  Stress 


Temperature 

°C 

MOR 

MPa 

WeibuJl 

Modulus 

Rate 

MPa/min. 

Number  of 

Samples 

30 

21 

787. 

6.5 

1870.0 

704 

730. 

6.0 

1.8 

5 

871 

818. 

8.5 

1.8 

5 

1038 

666. 

9.9 

1.5 

5 

1204 

485. 

8.0 

.8 

5 

1371 

313. 

11.2 

91.0 

15 

Reference  23 
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